Symmetry and localization of quantum walk induced by extra link in cycles 
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It is generally believed that the network structure has a profound impact on diverse dynamical 
processes taking place on networks. Tiny change in the structure may cause completely different 
dynamics. In this paper, we study the impact of single extra link on the coherent dynamics modeled 
by continuous-time quantum walks. For this purpose, we consider the continuous-time quantum walk 
on the cycle with an additional link. We find that the additional link in cycle indeed cause a very 
different dynamical behavior compared to the dynamical behavior on the cycle. We analytically 
treat this problem and calculate the Laplacian spectrum for the first time, and approximate the 
eigenvalues and eigenstates using the Chebyshev polynomial technique and perturbation theory. It 
is found that the probability evolution exhibits a similar behavior like the cycle if the exciton starts 
far away from the two ends of the added link. We explain this phenomenon by the eigenstate of 
the largest eigenvalue. We prove symmetry of the long-time averaged probabilities using the exact 
determinant equation for the eigenvalues expressed by Chebyshev polynomials. In addition, there 
is a significant localization when the exciton starts at one of the two ends of the extra link, we show 
that the localized probability is determined by the largest eigenvalue and there is a significant lower 
bound for it even in the limit of infinite system. Finally, we study the problem of trapping and 
show the survival probability also displays significant localization for some special values of network 
parameters, and we determine the conditions for the emergence of such localization. All our findings 
suggest that the different dynamics caused by the extra link in cycle is mainly determined by the 
largest eigenvalue and its corresponding eigenstate. We hope the Laplacian spectral analysis in this 
work provides a deeper understanding for the dynamics of quantum walks on networks. 
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I. INTRODUCTION 

The dynamic processes taking place in networks have 
attracted much attention in recent years [H-Q . It is gen- 
erally believed that network structure fundamentally in- 
fluences the dynamical processes on networks. Investiga- 
tion on such aspect could be done using the spectral anal- 
ysis and it has been shown that the dynamical behavior 
is related to the spectral properties of the networks [1, 0] . 
Examples include synchronization of coupled dynamical 
systems [sj, epidemic spreading percolation [tI], com- 
munity detection 8] , and others J5J . In several of these 
examples, the largest eigenvalue plays an important role 
in relevant dynamics. All these examples suggest that 
spectral property is crucial to understanding the dynam- 
ical processes taking place in networks. 

Quantum walks, as coherent dynamical process in net- 
work, have become a popular topic in the past few 
years The continuous interest in quantum walk 

can be attributed to its broad applications to many dis- 
tinct fields, such as polymer physics, solid state physics, 
biological physics, and quantum computation [ij, [l^ . 
In the literature there are two types of quan- 

tum walks: continuous-time and discrete-time quantum 
walks. It is shown that both type of quantum walks is 
closely related to the spectral properties of the Laplacian 
matrix of the network. Most of previous studies have 
studied quantum walks on some simple graphs, such as 



the line [1 
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, cycle [13, hvpercub e fisl . trees [l9l.[20|. 
I, ultrametric spaces |22| . threshold net- 
and other regular networks with simple topol- 
121. The quantum dynamics displays different 




behavior on different graphs, most of the conclusions 
hold solely in the particular geometry and how the struc- 
ture influences the dynamics is still unknown. Because 
quantum walks have potential applications in teleporta- 
tion and cryptography in the field of quantum compu- 
tation ^1^, it is clearly beneficial to investigate how the 
structure influences the dynamics of quantum walks. 

In this paper, we focus on continuous-time quantum 
walks (CTQWs) and study the impact of single extra link 
on its coherent dynamics. For this purpose, we consider 
the continuous-time quantum walk on the cycle with an 
additional link. The dynamics of continuous-time quan- 
tum walks on cycle is well known, the problem is analyti- 
cally solvable and directly related to quantum carpets in 
solid state physics The topology of cycle is highly 

symmetric and the quantum dynamics reflects such topo- 
logical symmetry. If one extra link is added into the cy- 
cle, the topological symmetry is broken and the structure 
has a small difference compared to the cycle. This en- 
ables us to treat the problem analytically and investigate 
how an additional link added to the network affects the 
dynamics. Similar study of the impact of single links 
on dynamical processes can be found in Refs. [24 l25j . 
where the authors study the impact of single link addi- 
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tion in percolation f23| and impact of single link failure 
in quantum walks j25|] respectively. In our case, in view 
of experimental implementation of CTQWs on cycles has 
been realized, we hope our study of the impact of single 
link (extra coupling) provides useful insight in quantum 
computation. 

The rest of the paper is organized as follows. Section II 
introduces the model of continuous-time quantum walks 
and defines the network structure. Section III investi- 
gates the eigenvalues and eigenstates of the Laplacian 
matrix (Hamiltonian). We obtain the determinant equa- 
tion for the eigenvalues expressed by Chebyshev polyno- 
mial, and calculate the largest eigenvalue on certain limit 
condition. The other eigenvalues and eigenstates are ob- 
tained by the perturbation theory. As we will show, like 
other dynamical processes, the largest eigenvalue and its 
eigenstate play an important role in the relevant dynam- 
ics. In Section IV, we study the characteristics of proba- 
bility evolution and distribution. The impact of the extra 
link depends on the initial exciton position, which can be 
well understood by the contribution of the eigenstate of 
the largest eigenvalue. The long time averaged proba- 
bilities shows significant symmetry and localization, and 
we give an explanation to these features using the deter- 
minant equation and largest eigenvalue. In Section V, 
we study the trapping process and the survival probabil- 
ities show significant localization for some special values 
of network parameters. We discuss this finding using the 
perturbation theory again. Conclusions and discussions 
are given in the last part. Sec. VI. 



II. CONTINUOUS-TIME QUATUM WALKS 
AND THE ADDITION OF LINK IN CYCLE 

A. Continuous-time quatum walks 

The coherent exciton transport on a connected net- 
work is modeled by the continuous-time quantum walks 
(CTQWs), which is obtained by replacing the classical 
transfer matrix by the Laplacian matrix , i.e., H = 
—T [m, The transfer matrix T relates to the Lapla- 
cian matrix by T = —7^, where for simplicity we assume 
the transmission rates 7 of all bonds to be equal and set 
7 = 1 in the following [ill, [13 • The Laplacian matrix 
A has nondiagonal elements Aij equal to —1 if nodes i 
and j are connected and otherwise. The diagonal ele- 
ments An equal to degree of node i, i.e., An = ki. The 
states |j) endowed with the node j of the network form 
a complete, ortho-normalised basis set, which span the 
whole accessible Hilbert space. The time evolution of a 
state |j) starting at time to is given by \j, t) = U {t, io) I j) , 
where U(t,to) — exp[—iH{t — to)] is the quantum me- 
chanical time evolution operator. The transition ampli- 
tude ak,j{t) from state \ j) at time to state \k) at time 
t reads ak.i(t) = {k\U{t,0)\j) and obeys Schrodingers 
equation |lll . |l7[. The classical and quantum transi- 
tion probabilities to go from the state \j) at time to 



the state |fc) at time t are given by Pk At) = {k\e 
and TTkAt) = WAt)? = l(fc|e-**^|j)|' [III El, respec- 
tively. Using En and jvPn) to represent the nth eigenvalue 
and orthonormalized eigenstate of H, the quantum am- 
plitudes between two nodes can be written as [lli] 



n.At) = \<^k,{t)\' ^\Y.e-^'^-{k\^n){^n\j)\' 

n 

^J2e-^'^^--^'\k\^n){^n\j)U\^i){^i\k). 

n,l 

For finite networks, TT^jit) do not decay ad infinitum 
but at some time fluctuates about a constant value. This 
value is determined by the long time average of tt^ j (t) 

1 

Xkj = lim - / TTk,j{t)dt 

= Y,{k\^n){^n\j)mi){^l\k) 

, ,T (3) 
X lim - / e-**(^"-^')dt 

= Y,S{En,Ei){k\^n) Ij) (j I I fc) . 

where S{En,Ei) takes value 1 if En equals to Ei and 
otherwise. To calculate the exact analytical expressions 
for TTk.j (t) and Xk,j , all the eigenvalues En and eigenstates 
I'^n) of the Laplacian matrix are required. If all the 
eigenvalues {Ei \i = 1,2,..., A''} are distinct, i.e., all the 
eigenvalues are not degenerated, Eq. ([3]) can be simplified 
as, 

Xfc..=ElW*")M(*"|j')l'- (4) 

n 

The eigenvalues for the cycle graph is twofold degener- 
ated. However, as we will show, the addition of single 
link causes the degeneracy to disappear. Therefore we 
can use Eq. ([4]) to calculate the probability distribution 
for our model. The transition prababilities (See Eqs. 

are closely related to En and |vl/n), it is crucial to 
calculate the Laplacian eigenspectrum and we will do this 
in Sec. III. 



B. Addition of link in cycle 

Now we specify the topology of the cycle with an extra 
link. First we construct a cycle of size N where each node 
connected to its two nearest-neighbor nodes, then we con- 
nect two node of certain distance with an additional link. 
The topology is completely determined by the distance 
and size of the cycle. For the sake of simplicity, we assign 
a consecutive number from 1 to TV for each neighbored 
node in the cycle. The additional link connects node 1 
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the second largest eigenvalue does not converge to zero 
when the size of the system goes to infinity, whereas the 
other nearest eigenvalue gaps tend to zero in the limit of 
infinite system. This suggests the eigenvalue spectra is 
continuous except the largest eigenvalue. Such charac- 
teristic feature of eigenvalue spectra is useful for under- 
standing the coherent dynamical behavior. As we will 
show, the largest eigenvalue plays an important role in 
the dynamics and gives significant contribution to the 
localized probabilities. Hence in the following, we will 
try to obtain Laplacian eigenspectrum and determine the 
largest eigenvalue and its corresponding eigenstate. 



A. Determinant equation for the eigenvalues 



FIG. 1. Topology of the cycle with an extra link. The addi- 
tional link connects node 1 and m = 3 in the cycle. According 
to the symmetry of topology, the graph can be projected into 
ID chain of size N' ^ \N/2] (the notation \N/2] denotes the 
integer part of N/2). 



and node m in the cycle. Thus the structure, denoted 
by G{N,'m), are characterized by the network size N 
and connecting node m {m G [3, iV — 1]). Here, we use 
d{l,m) to denote the (shortest) distance between node 
I and node m in the cycle (without an extra link), the 
graph can also be characterized by G{N,d{l,m)). The 
structure of G{N, m — 3) is illustrated in Fig. [TJ 

It is interesting to note that the network considered 
also has a symmetric structure, and the symmetry axis 
lies at the central position between node 1 and m. The 
structure can be mapped into a ID chain of size N' « 
(the notation denotes the integer part of 

N/2) if we make a horizontal projection of the graph 
(See Fig. [T]). Such projection sheds some light on the 
implicit relationship between the two structures. For ex- 
ample, as we show in Appendix C, the eigenvalues of the 
Laplacian matrix of the projected ID chain also belong 
to eigenvalues of our model. 



III. LAPLACIAN EIGENVALUES AND 
EIGENSTATES 

Since the transition probabilities are determined by the 
Laplacian eigenvalues and eigenstates (see Eqs. (I2])-(|4])), 
a detailed analysis to the eigenvalues and eigenstates will 
be helpful for the problem. In this section, we will study 
the Laplacian eigenvalues and eigenstates in detail. 

Figure [Ha) shows the eigenvalues obtained by numeri- 
cal diagonalizing the Laplacian matrix using the software 
Mathematica 7.0. The eigenvalues are ranked in ascend- 
ing order for networks of size = 100 with m = 3, 
m = 5 and m = 10. It is interesting to note that the 
largest eigenvalue is larger than 4 and isolated from the 
other eigenvalues (the other eigenvalues are less than 4) , 
this means the gap between the largest eigenvalue and 



We start our analysis on the eigen equation of the 
Hamilton(Laplacian matrix). The Laplacian matrix H 
of G{N, m) (rn ^ 3) takes the following form: 



m\j) 



3, if i = J = 1, or i — j ^ m 

2, a i ~ j ^ 1, and i =^ j ^ m 

— 1, if i and j connected 

0, otherwise. 



(5) 

According to the eigen equation H\^) = E\'^), suppose 
the eigenstate 1^") can be expressed as 



N 

1=1 



(6) 



The eigen equation can be decomposed into the following 
N linear equations, 

iXi - X2 - Xm - ^ Exi, (7) 

—Xj^i + 2xj — Xj+i — Exj, 1 < j < m (8) 

Xl Xm—i -\- 3x^yi — EXjji^ (9) 

^Xj-i + 2xj — Xj+i = Exj, m < j < N (10) 

— Xn-~1+2Xn^Xi—EXn. (11) 

Eq. (ISl) can be rewritten as (2 — E)xj = xj-i + Xj+i. 
This is similar to the recursive definition of the Cheby- 
shev polynomials of the second kind (See Appendix A). 
Noting the recursive relations and setting {2 — E) = 2x 
in the definition of Chebyshev polynomials, the variables 
xi,X2, Xm-2 in Eq. can be expressed as a function 
of Xm-i and Xm, 

Xj-l = Um-j{x)Xrn-l-Um-i-j{x)Xm- 1 < j < m (12) 

Analogously, Eqs. (|10p and dTTI) can be written as a func- 
tion of a;]v and xi using the Chebyshev polynomials, 

Xj-i — U„+i-j{x)xj^ — UN-j{x)xi. m < j ^ N (13) 

Substitute E ^ 2 - 2x into Eqs. dZ]) and (O, we obtain, 

3xi — X2 — Xm — Xjv = (2 — 2x)xi, (14) 
^{2-2x)x^. (15) 
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In Eq. ([12]), Xi = Um-2{x)Xm-l - Um-3{x)Xm, X2 = 

Um,-3ix)xm-i — Um-4{x)xm. Replacing the variables Xi 
and X2 in Eqs. (fT4|) and (fT5|) . we get, 



[{2X + l)Urn-2{x) - U„,-3{x)] 

[Um-i{x) - (2x + l)U„^-z{x) - 1] 

Xyn+1 = -(1 + ?/m-2(a;))x™_i + 
(2x + 1 + [/™_3(a;))xTO. 



(16) 



(17) 



Substituting Xi = Um-2{x)xm-i ~ Um-3{x)xm and Xm 
in Eq. (fT6| into Eq. (fT3|) for j = rn + 1 and j = m + 2, 
we obtain. 



CiXm-l + C2X„ 



0, 



(18) 



where ci = C/jv-m(a;)[(2a; + l)[/,„_2(a;) - Uyn-z{x)] ~ 
Ur,-m-iix)U„i-2{x) and C2 = [/jv-m(a;)[[/„_4(a;)- (2a; + 
l)t/m-3(a;) - 1] + UK-,n-i{x)U,n-3{x) - 1, and 

Xm+i = {C/jv-m-i(a;)[(2x + l)t/„_2(a;) - Um-3{x)] 

- UN-m-2{x)Um-2ix)}Xm-l + 

{Uj,-rn-l{x)[Um-i{x) - (2a; + l)Um-3ix) - 1] 
+ UN-ni-2{x)U,n-3{x)}Xm- 

(19) 

Combine Eq. (fTT)) and Eq. (|T9l) , we get another equation 

for Xm^i and Xm, 



C3Xm-l + C4X„ 



0, 



(20) 



where C3 = [/jv_m_i(a;)[(2a; + l)?7TO_2(a;) - 

Um-3{x)] - Ur,-m-2{x)Um-2{x) + Um-2{x) + 1 and 
C4 = UM-m-l{x)[Um-4.{x) ~ {2x + l)Um-3{x) ^ 1] + 

Ur,-m~2ix)Um~3{x) - Uyn-3{x) - {2x + 1). Thus we 
have got two equations for Xm-i and Xm, Eq. (fT8|) 
and Eq. (piU) . The two equations should have nonzero 
solutions, which leads to. 



C1C4 — C2C3 = 0. 



(21) 



In the Appendix B, we show the substraction of the 
products can be simplified to be a much simple form in 
Eq. (jB7[) . thus we have obtained determinant equation 
for the eigenvalues, 

1 + U^^,n{x) + U,n^2{x) - [/„-l(x) - T^{x) = 0, (22) 

where T^{x) and Um-i{x) arc Chebyshev polynomials of 
the first kind and the second kind, respectively. Solving 
the above equation, we can get all the Laplacian eigen- 
values. This equation is useful to determine the largest 
eigenvalue and interpret the symmetric structure of the 
transition probabilities. 



B. The largest eigenvalue 

The largest eigenvalue can be determined using 
Eq. (1221) under certain limit conditions. As we have 



shown, the largest eigenvalue is larger than 4, this cor- 
responds to the solution a; < — 1 in Eq. ([2^ . Note that 
the Chebyshev polynomials is divergent for |a;| > 1 in the 
limit of infinite order, Eq. (P^ divided by Tm{x) leads 
to. 



Ui^-yn{x) ^ Uyn-2{x) _ U„^i{x) 

Tn{x) Tm{x) Tm{x) 



1 = 0. 



(23) 



1 less 
'V2. Thus 



For large size of system and long range coupling, i.e., 
N is large and the distance d{l,m) between node 1 and 
node m is not small, if we apply the asymptotic solution 
of the Chebyshev polynomials (See Eq. (jA5l) in Appendix 
A), the first two terms of the above equation equal to 0, 
_Un,M^ _ 2 — ^ ^ ^i^i^j^ igj^jg to zo = -1 - V2 

(the other solutions do not satisfy z = x 

than —1 when x < —1) and xq = "22^ 
the largest eigenvalue equals to a constant value E™''^^ = 
2 — 2xo = 2 + 2\/2 under this limit condition. 

To test the above prediction, we plot the largest eigen- 
value as a function of N and m in Fig. [UJb). As we 
can see, the largest eigenvalue converges to a constant 
value as iV or TO increases. For moderate range coupling, 
e.g., TO — 5, the largest eigenvalue is close to the ana- 
lytical prediction (2 -I- 2^/2). The constant value of the 
largest eigenvalue suggests that its corresponding eigen- 
state (eigenvector) approaches to certain stationary dis- 
tribution. Here, the obtained largest eigenvalue is useful 
to consider its corresponding eigenstate, which we will 
show in the following. 



C. Eigenstate of the largest eigenvalue 

The components of eigenstate corresponding to the 
largest eigenvalue are shown in Fig. [IJc). The straight 
line in the linear-log plot indicates that the components 
display an exponential decay, i.e., \x'j\ ~ Zq'' . The dis- 
tribution of the components of eigenstate is symmetric. 



V j e [1,to]; 



I, VjG 



[m + 1,N]. As we will prove in Sec. IV (subsection B) 
and Appendix C, this symmetric behavior is true for all 
the eigenstates. 

To obtain an approximate formula for the eigenstate of 
the largest eigenvalue, we use Eq. (fT8|) to present Xm~i as 
a function of Xm, and substitute it into Eq. (|16|) to write 
X as a function of Xm ■ Using the Chebyshev polynomial 
identity (jA7[) in Appendix A, Eq. (IT6l) can be written as 

Xn = [Ujn-lix) + Um-2{x)]x„i-l-[Um-2ix) + Um-3ix) + 

l]xm- According to Eq. (fT8| . x„i-i = —C2Xm/ci and 
noting the simplified ci , C2 in (jBl[) and (IB2I) in Appendix 

A, we obtain x^ = -f-Xm (See the detailed derivation in 
Eq. (jBSP in Appendix B). Substituting Xm-i and x^ into 
Eqs. (fT^ and (fT5|) and noting xi = ~x,n for the largest 
eigenvalue (See Eq. (IC12I) in Appendix C), Eqs. (fT^ and 
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FIG. 2. (Color online) (a)Eigenvalues Ej, arranged in ascending order, for the networks of TV = 100 with m = 3, m = 5 and 
m — 10. The largest eigenvalue is larger than 4 and isolated from the other eigenvalues, (bl) shows the largest eigenvalue E™'"'^ 
as a function of system size A*' for m = 3, m = 5 and m = 10. (b2) shows the largest eigenvalue E™'°'^ as a function of m for 
fixed network size A'^ = ICQ. In large system and long range coupling, i.e., the distance between node m and node 1 (d(l,m)) 
is not small (the plat region in the plot), the largest eigenvalue approaches to constant value (2 + 2\/2). (c) Components of 
eigenstate of the largest eigenvalue \x'j\ as a function of j for networks of size A'' = 100 with m = 20 (red dots) and m = 50 
(black squares). The decay is exponential and there are two maximal points (ja;'il and la::^!), the distribution is symmetric, 
\x'j\ = |a;^+i-j|, V j G [l,m], \x'j\ = \x'ra+N+i-j\, V j € [m + 1, A^]. 



(fT3|l can be written as, 

,C2 



Ura-]-l{x) + Um-2-j{:x)]Xjn,i £ [l,77l] 



Cl 



Xj = [—UN-jix) + UN-j~iix)]xm- j G [m,N] 

^1 



(24) 



For the largest eigenvalue, x = xa = — V2, the ra- 
tios C2/C1 and C2/C1 approach to the constant value 
—Zq^ = (v^ — 1) for large system and long range cou- 
phng. Therefore, Eq. (IM)) can be recasted as, 

^'j ~ [^0^ • Um-l-jixo) - U,n-2-j{xQ)]x[„ 



J — m f 



TO+I 



c^- « [-Zq ^ • UM-j{xo) + U,^-i-.j{xo)]x'„ 



(25) 







Xm- .1 e 



.{m + N + l) 



^,N] 



The above formula only describes the right part of 
the symmetric distribution of the eigenstate (See the 
points on the dashed lines in Fig. (He)). The decay 
exponents in Eq. p5|) only depend on the (shortest) 
distance between the node pairs (j, 1) and (j, m), i.e., 
dj = miii{d{j,l),d{j,m)}. Eq. ([^5]) can be summarized 
in a much simple form using the shortest distance dj, 

x' w ±z^'^'x'^^, dj = min{d(j, l),d{j, m)) (26) 



In the calculation, the eigenstate should be normalized. 



N 



r-l |2 



1. For large system and long-range couple, 



the summation has a week dependence on N and m, thus 
the summation of finite geometric series approaches to 
the summation of the infinite geometric series, i.e., 

N m N 

j=l j = l j=m+l 



'')\xL\' 



(27) 



^(2E-o'"''+2E-o 
--2V2\x'„f = 1, 



which leads to \x[\ = \x'„^\ — 2^'^/'*. For large sys- 
tems and long range coupling, the distribution of the 
eigenstate of the largest eigenvalue approaches to a sta- 
tionary exponential decay with the maximal component 
\x[\ = \x'^\ = 2-3/4. Fig.EJc), we plot the analytical 
predictions in Eq. (|26p , which agrees well with the exact 
numerical result. 

We would like to point out that the exponential distri- 
bution of the components of the eigenstate is an interest- 
ing characteristic for the largest eigenvalue. For the other 
eigenstates, the components do not display this feature. 
This could be interpreted by the striking different behav- 
ior of the Chebyshev polynomials for |a;| > 1 and |x| < 1. 
For I a; I ^ 1, the Chebyshev polynomials show finite reg- 
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ular oscillations, but for \x\ > 1, Chebyshev polynomials 
show exponential growth. The largest eigenvalue and its 
eigenstate play a significant role in the coherent dynam- 
ics, we will discuss this in the next section. 



D. Other eigenvalues and eigenstates 

We have determined the largest eigenvalue and its cor- 
responding eigenstate using the determinant equation 
()22|) . However, it is not intuitionistic to determine the 
other eigenvalues and eigenstates by Eq. ([22l) . Here, we 
use the perturbation theory to get the approximate re- 
sults for the other eigenvalues and eigenstates. 

According to perturbation theory, the Hamilton 
(Laplacian matrix) of the cycle with an additional link H 
can be divided into two parts: the unperturbed Hamil- 
ton and perturbed Hamilton H' {H = + H' , 
H\-^) = £'!*)). The unperturbed Hamilton is the 
Laplacian matrix for the cycle, whose eigenvalues and 



eigenstates are well known: H'-^'>\^jj 
= 2-2cos0„ , ivi"^) = 



[- rf 1 ' rf 1]' ^« = 2mT/N). The perturbed Hamilton H' 
can be written as H' — |1)(1| -I- \m){m\ ~ |l)(m| — |m)(l|. 

The eigenvalues of the cycle E^i ^ are twofold degenerated 
except the minimal eigenvalue and maximal eigenvalue 
4. The first order corrections for eigenvalues and 4 



are Ej^'' = and E\^>^ = 2/iV[l + (-1)™] respectively 
(maximal eigenvalue 4 only exists for even N). The zero 
order approximation of the eigenstates corresponding to 
and 4 are equal to the original unperturbed eigenstates 

of z.e., 1*^°)) = iV^r) - yVNEU 1*^/2) = 

l^/J'/2) = V\/^E^i(-l)^|j)- For the twofold degener- 
ated states, the degeneracy disappears when an extra link 
is added. To approximate these eigenvalues and eigen- 
states, we apply degenerate perturbation theory to cal- 
culate the first order corrections of the eigenvalues and 
zero order approximation of the eigenstates. 

We diagonalize the the perturbed Hamilton H' us- 
ing the unperturbed degenerated eigenstates IV'n'^) and 
lij-n)- The perturbed Hamilton H' can be written as. 



The corresponding eigenstates for En^ and isi"'^ are, 



V2\ 



i{m+l)e, 
1 



1 



„i(m+l)e„ 

1 



(30) 



The zero order approximation sof the eigenstates are 
given by 



|*W)=if„(l)IV'i°^)+i^n(2)|Vi°I) 



I* 



(0)n 



=/^_„(l)|VA("))+i^-„(2)|^L"^) 



(31) 



/,(0)\ 



Substitute the Bloch states iV'ri / — 
into the above equation, we obtain. 



1 



I* 



(0)n 



N 
N 

=v 



1 



(e 



,-iO-m-l)fl 



-i{j-m-l)( 



Ob), 



(32) 



Thus we have obtained the first order corrections of the 
degenerated eigenvalues and the zero order approxima- 
tion of the corresponding eigenstates(See Eqs. and 
([5^ 1. Using these eigenvalues and eigenstates, we can 
calculate the transition probabilities in Eqs. ([2])- ([4]). 

It is worth mentioning that the perturbation theory 
can not be used to calculate the largest eigenvalue, since 
the eigenvalue gaps do not converge to zero in this case. 
However, the analytical predicted eigenvalues by pertur- 
bation theory (See Eq. (PU)) ) can also be obtained by 
Taylor series expansion of the Chebyshev polynomials. 
In Appendix D, we approximate the other eigenvalues 
using the determinant equation (1^^ . which are consis- 
tent with the predictions in Eq. (P^ . 



IV. PROBABILITY EVOLUTION AND 
DISTRIBUTION 

In this section, we consider the probability evolution 
and probability distribution. We will explain the ob- 
served phenomenon using the obtained eigenvalues and 
eigenstates. 



H' 



where A 




[l,rfl - A] 



(28) 



if N is odd 



where i?,;„ = H\_„ = {iP^r?^\H'\^i°^) = f [l-cos(m- 

e™^"]. The first order corrections of the eigenvalues are 
given by the two eigenvalues in Eq. 



Ei'^=0, E 



(1) 



N 



[1 — cos(to — 1)6„ 



(29) 



A. Probability Evolution 

The time evolution of the probability is given by 
Eq. ([2]). Here, we focus on the return probability TTjj{t) 
(set k — j m Eq. ([2])). The decay behavior is used to 



quantify the efficiency of the transport [11| . In our case 
the return probability ttjj (t) depends on the starting po- 
sition j. We study the time evolution of the return prob- 
ability T^j,j{t) on different starting position j. The be- 
havior of TTjj (t) versus t are plotted in Fig. |31 It is found 
that TTjj{t) shows different time range for scaling behav- 
ior Trjj{t) ~ t~^: For the starting positions far from the 
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plitude comes from the other eigenvalues, which leads to 
the scaling behavior resembling to the case in cycles. If 
the exciton starts at the center nodes (nodes at symme- 
try axis) j = {771 + l)/2 (m is odd) or j = (iV + m + l)/2 
{N + 771 is odd), the corresponding components of the 
largest eigenvalue's eigenstate are (See Eqs. (|C13p - 
(IC14I) in Appendix C). In this particular case, the largest 
eigenvalue gives no contribution to the amplitude, the 
probabilities are exact the same as in the cycles, thus the 
additional link has no impact on the dynamics. If the ex- 
citon starts at the positions near the connected nodes 1 
and m, the components of the largest eigenvalue's eigen- 
state become larger and give contribution to the ampli- 
tude, thus the scaling behavior disappears gradually (See 
Fig. |31[b) (c)). If the exciton starts at nodes 1 or m, the 
components of the eigenstate of the largest eigenstate (xi 
or Xm) give large contribution to the amplitude, the time 
range of the scaling behavior is very short (See Fig.[3l^d)). 
In this case, the return probabilities do not decay ad in- 
finitum but fluctuate about a constant value, suggesting 
significant localization induced by the largest eigenvalue. 
We will study this feature in detail in the following. 



B. Probability Distribution 



FIG. 3. Return probabilities 7rjj{t) on a netwok of A*' = 
100 and m = 21 for different starting positions j. (a)The 
exciton starts at the center node (symmetry axis) between 1 
and m, i.e., j = {m + l)/2 = 11 or j = (iV -I- m -I- l)/2 = 61. 
(b)Starting at j — 41. (c)Starting at j = 31. (d)Starting 
at J = 1 or J = m = 21. The time range of scaling 
becomes short when the starting position near the two ends 
of the added link. 



connected nodes 1 and rn, the scaling range of time is 
much larger than that of the walk starting near the nodes 
1 and 717, (Compare scaling in the four figures). It is well 
known that for the cycle, TTjj{t) ~ ^o(2t) ^ t^^. This 
suggests that the additional link gives small influence on 
the dynamics when the walk starts at outlying positions, 
and gives larger influence for exciton near nodes 1 and m. 
If the walk starts at the center nodes (nodes at symmetry 
axis) j = {7n + l)/2 (m is odd) or j = {N + 77i + l)/2 
{N + 771 is odd) (See Fig. [3] (a)), the probabilities are 
exact the same as the probabilities in cycle of the same 
size, in this case the additional link has no impact on the 
dynamics. 

The observed phenomenon can be well understood us- 
ing the eigenstates obtained in the previous section. If 
the exciton starts at positions far away from 1 and m, the 
components of the eigenstate of the largest eigenvalue are 
very small due to the exponential decay of the compo- 
nents. Thus the contribution from the largest eigenvalue 
in Eq. ([T|) is very small, the main contribution to the am- 



In this section, we consider the long time averaged 
probabilities. We find that these limiting probabilities 
display a symmetric structure, symmetric nodes k and I 
{k + l = 7n + loTk + l = N + 7n + l, see the two nodes 
in the horizontal projection in Fig.[T]) have the same lim- 
iting probabilities. This suggests that the probability of 
finding the exciton at a certain node k is equal to the 
probability of finding the exciton at the symmetric po- 
sition / = (rn + 1) — fc (or / = N + 7n + 1 — k), i.e., 
Xk,j = Xm+i-fcj V j,k (or Xfcj = XN+m+i-k,j V j,k). 
This symmetric structure can be understood as a result 
of the axis symmetry of graph. Here, we give a rigorous 
mathematical proof for this symmetric structure using 
Eq. Q. In Eq. (U), the limiting probabilities are only 
dependent on the eigenstates. Suppose the nth eigen- 
state is expanded as |^„) = J2jLi ^'j^^lj)^ ttie probabili- 
ties in Eq. ^ can be written as Xk,j = J2n l-^j"^ l^l^l"^ 1^- 
In Appendix C, we prove that the components of eigen- 
state are symmetric, \xj\ = |a:„i+i_j|, V j € [1,to]; 
\xj\ = \Xm+N+i-j\, V j e [m + 1, A^]. Here, no super- 
script for Xj means the identity holds for all the eigen- 
states (V n) . The symmetric structure of the eigenstates 
leads to the symmetric limiting probabilities. This is one 
of the main conclusions in our paper. 

Now we try to use the eigenstates to find analytical ap- 
proximate results for the limiting probabilities. The con- 
tribution in Eq. ((4]) comes from three parts: the eigen- 
state of the largest eigenvalue, non-degenerated eigen- 
states d^o"^) and |5'^/2)) ^^'^ degenerated eigenstates. 



8 



0.03 




20 40 60 80 100 



0.02 



0.00 



0.15 







1,1,1 


. , .(c) 



20 40 k 60 80 100 




20 40 k 60 80 100 



FIG. 4. (Color online) Long-time averaged probability distri- 
bution Xk,j on a network of TV = 100 and m = 21 for different 
starting positions j (compare Fig.[3l). (a)The exciton starts at 
the center node between 1 and m, i.e., j — (m-|-l)/2 = 11 and 
j = (A'' + m + l)/2 = 61. (b)Starting at j = 41. (c)Starting 
at j = 31. (d)Starting at j = 1 and j = m = 21. The black 
squares are the exact numerical results obtained by diago- 
nalizing the Laplacian matrix. The dashed red lines are the 
analytical predictions in Eq. ((34)|. 



Therefore Eq. (jlj transforms into 

n=l 
[Af/21-A 



(33) 



where A takes value 1 for even N and for odd N. Sub- 
stituting the degenerated eigenstates of Eq. (l32t into the 
above equation and noting the trigonometric function for- 
mula J2Ui cos ^ = i (sin (22^^ / sin ^ - 1) , Eq. §^ 
can be simplified as, 



Xk. 



1 f sin27r(j-fc)(l- A) 



sm 



N 



(34) 



sin27r(j + k - m - 1)(1 - ^) 

o.in 27r(j+fc-m-l) 
bin ^ 



Fig. |4] shows the probability distribution for different 
starting positions. It is evident that the exact numerical 
results agree well with analytical predictions in Eq. ([M]) . 
For exciton starting sd j = 1 and j — m, there is a 
high probability to find the exciton at j = 1 and j — 
TO, suggesting significant localization in the dynamics. 
For excitons starting at center node j = (to -I- l)/2 (to 
is odd) and j = {N + m + l)/2 {N + to is odd), the 
probability distribution is exact the same as the cycles. 
This could be explained by Eq. ([M)) as follows: For j = 
(to + l)/2 or j = {N + m + l)/2, \x'j\ = (See Eqs. 
(|C13l) - (jC14p in Appendix C), the contribution from the 
largest eigenstate's eigenstate is (See the first term in 
Eq. ([M]) ). the localized probability Xkj mainly depends 
on the trigonometric function in Eq. (j34p . For return 
probabilities with j — k = (to -|- 1)/2 (to is odd) or 
j = k = {N + m + l)/2 {N + m is odd), limits of the 
trigonometric function ratio in the big bracket equal to 
N — X, thus Eq. (IMl) becomes as. 



Xkj 



(2jV-l-A) 



J 



m+1 (N+m+l) 
2 ' 2 



(N-l-X) . / , . , _ rn+l (jV+m+1) 

Ar2 J J -/- i^T 1^ — 2' 2 



which are consistent with the results in Ref. [2f 



(35) 
For this 



particular case, j, fc = { ^22±i ^ L^Y±E±il } ^ the probability 
distribution is exact the same as the cycles, this indicates 
the adding of link has no impact on the dynamics for 
center node excitations {k or j at the symmetry axis. See 
Fig. [1]). Here we provide a theoretical interpretation for 
this phenomena. However, if to and N + m are not odd 
numbers, this particular case does not exist. 

As we have shown, the return probabilities Xj,j h^-s a 
high value at j = 1 and j = to. Such a strong localization 
may vanish as the size of the system becomes large. In 
order to investigate the dependence of localized probabil- 
ity Xi,i (or Xm,m) on the size of the system, we plot xi,i 
(or Xm,m) as a function of N for different values of to in 
Fig. [5l We find that the return probability converges to 
a constant value when N tends to infinity. This constant 
value can be calculated using Eq. ([M)) . For infinite sys- 
tems, xi,i (or Xm,m) equals to the first term of Eq. ([M)) 
(the other terms are in the limit of — oo). There- 
fore, for infinite systems, the limiting probabilities are 
determined by the components of the eigenstate of the 
largest eigenvalue. Noting that \x'i\ — — 2"^^/^, we 
get 



Xi,i = Xi, 



Xm.l Xji 



„' |4 



(36) 



Therefore xi,i (or Xm,m) converges to the 1/8, indicating 
there is a significant lower bound for the return prob- 
abilities. This result has not been previously reported 
and suggests that quantum walks display nontrivial lo- 
calization on nonsymmetric system. Since all the com- 
ponents of the largest eigenvalue's eigenstate are known 
(See Eqs. ((25ll and p6|). the limiting probabilities Xkj 



have a lower bound given by \x' 
sion of a;' in Eq. (f^ . 



. Using the expres- 

|^g|-2(<i,+d.)|^^ 
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authors consider a system of N nodes and among them 
M are traps (M < N). The trapped nodes are denoted 
them by to', so that m' G Al. The new Hamiltonian 
of the system is H — Hq + iF, where Hq is the orig- 
inal Hamiltonian without traps and iV is the trapping 
operator. T has m purely imaginary diagonal elements 
Trn'm' at the trap nodes and assumed to be equal for all 
TO {Tm'Tji' = r > 0). See Ref. [Ill for details. The new 
Hamiltonian is non-hermitian and has N complex eigen- 
values and eigenstates {Ei, |^/)} {I — 1,2,..., N). Then 
the quantum transition probability is 



20 40 N 60 80 100 



it)^\a,At)\' = \Y.' 



itE, 



{k\^i){^m\ (37) 



FIG. 5. (Color online) Return probabilities and x™,m ^ 
a function of the network size A'' for m = 3 (black squares), 
m — 5 (red dots) and m = 10 (blue triangles). All the return 
probabilities approach to the predicted constant value 1/8 
(The horizontal dotted line). 



^\zq\^^'^''-^^'^''\ The lower bound for the limiting proba- 
bility is determined by the shortest distance dj and d^. 
For instance, for j = 1 and k = 2, dj = 0, dk — 1, 
> s{i+V2y ' J = ™ and fc = iV - 1, dj = 0, 

dk 



2, Xk,, > 



etc. This dynamical fea- 



ture is completely different from the case of the cycles. 
For continuous-time quantum walks on a ID cycles of 
N nodes, all the limiting probabilities Xk,j are when 
— ^ oo, due to all the components of the eigenstates 
approaches to in this case. So there is no localization 
on infinite regular cycles. For our case, the localization is 
induced by the largest eigenvalue and size of the system 
N, which differs from the ID regular cycles where the 
localization is only induced by the size of the system. To 
this end, the localization in our case is much more es- 
sential and may have potential applications in quantum 
information science. 



V. TRAPPING 

An im por tant process related to random walk is trap- 
ping [27l. l28l| . Trapping problems have been widely stud- 
ied in the frame of physical chemistry, as part of the gen- 
eral reaction-diffusion scheme poj . Previous work has 
been devoted to the trapping problem on discrete-time 
random walks (sol . [3l| . However, even in its simplest 
form, trapping was shown to yield a rich diversity of re- 
sults, with varying behavior over different geometries, di- 
mension, and time regimes \3^. The main physical quan- 
tity related to trapping process is the survival probabil- 
ity, which denotes the probability that a particle survives 
during the walk in a space with traps. 

In this paper, we consider trapping using the approach 
based on time dependent perturbation theory and adopt 
the methodology proposed in Ref. [11]. In Ref. 32 1, the 



where {I — 1,2,...,N) is the conjugate eigenstates 
of the new Hamiltonian. Equation ([57| depends on the 
initially excited node j. The average survival probability 
over all initial nodes j and all final nodes k, neither of 
them being a trap node, is given by [32| . 



UMit) 



1 



N - M 



(38) 



For intermediate and long times and a small number of 
trap nodes, H^,/(t) is mainly a sum of exponentially de- 
caying terms ^32i] : 



JV 



N~ M 



^exp(-27,i). 



(39) 



1=1 



where 7/ is the imaginary part of the eigenvalue Ei . 

For the cycle with an additional link, we focus on the 
case that only one trap exists in the graph and the trap is 
located on node 1. Fig. [6] shows the survival probability 
HM(i) vs i for AT = 100 (a) and N = 101 (b). We find 
that llM{t) decays very slowly and converges to a cer- 
tain constant value (nonzero) for some special network 
parameters N and m. Meanwhile, for the other values of 
N and m, IlM{t) decays slowly and converges to 0. For 
the network of size N = 100 (See Fig.lHlIa)), HmW con- 
verges to nonzero constants for to = 3, to = 6, to = 11, 
TO = 26 and to = 51, whereas IlM{t) approaches to 
for other values of m. For the network of size = 101 
(See Fig. Efb)), IlM{t) approaches to for all values of 
TO. The nonzero threshold of the survival probability 
suggests that the trapping process displays significant lo- 
calization for some special values of network parameters. 
This result is interesting since in previous studies the sur- 
vival probability converges to [llj, [s^, indicating the 
excitation eventually returns to the traps. The nonzero 
threshold of the survival probability implies that the ex- 
citation does not eventually returns to the traps [33| , and 
the survival probability is smaller than 1/2. 

To give an explanation to the localized survival prob- 
ability, we analysis the Laplacian eigenvalues using the 
perturbation theory. In Eq. ([55]) . the behavior of Hjv,/(t) 
is determined by the the imaginary part of the eigenvalue 
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FIG. 6. Survival probabilities nA/(t) for = 100 (a) and 
= 101 (b) for different values of m. In the calculation, 
the trap is localized at node 1 and F = 1. We find that if 
2(m-i) integral, IlM{t) converges to (m — 1 — ^)/{N — 1). 
Otherwise, nji/(t) approaches to 0. 



(ji). If 7i is nonzero (7; > 0), exp(— 27;^) converges to 
at long time scale. On the contrary, if ji is 0, exp(— 27;t) 
equals to 1. The main contribution for the localized sur- 
vival probability comes from the real eigenvalues (7; = 0), 
thus the lower bound of Ili\i{t) depends on the number 
of real eigenvalues (7; — 0). The eigenvalues and eigen- 
states for cycle with an additional link has been analyzed 
in Sec. Ill, the trapping operator can be regarded as the 
perturbed operator. According to perturbation theory, 
the first order corrections of the eigenvalues of the non- 
hermitian Hamiltonian are given by 



i?w = (f(°)|zr|*(")) = zr|(i|M^(''))|=^ 



(40) 



Utilizing the zero order approximation of the eigenstates 
in Eq. (|32|) . the first order correction of the eigenvalue 
can be written as En^ = ir-^[l ± cos(to — 1)6'„]. 7; = 



requires the first order correction E., 

2n(m— 1)71 



(1) 



equals to 0, thus 



(m - l)9n = ilH!^ = fcTT (n e [1, \N/2] - A], fc > 0). 
The number of integer solutions of this equation gives 
number of real eigenvalues 7/ — 0. If 2{m-i) hrtegral, 
the number of integer solutions is (m— 1 — A). Therefore, 
the lower bound of Hm {t) equals to (m — 1 ~ X)/ {N — 1) 
(< i). If N is odd, no integer solutions exist, thus njv/(t) 
converges to for odd size of networks. Thus we succeed 
in explaining the localized survival probability. 



VI. CONCLUSIONS AND DISCUSSIONS 

In summary, we consider the continuous-time quantum 
walk on the cycle with an additional link. We analytically 



treat this problem and and approximate the Laplacian 
eigenvalues and eigenstates by the Chebyshev polynomial 
technique and perturbation theory for the first time. We 
find that the probability evolution exhibits a similar be- 
havior like the cycle if the exciton starts far away from 
the two ends of the added link. The distribution of the 
long-time limiting probabilities display symmetric struc- 
ture, we prove this symmetry using the exact determi- 
nant equation for the eigenvalues expressed by Cheby- 
shev polynomials. In addition, the quantum dynamics 
exhibit significant localization when the walk starts at 
the two ends of the extra link, we show that the localized 
probability is determined by the largest Laplacian eigen- 
value and there is a significant lower bound for it even in 
the limit of infinite system. Finally, we study the prob- 
lem of trapping and show the survival probability also 
displays significant localization for some special values of 
network parameters, we succeed in explaining this phe- 
nomenon and determining the conditions for the emer- 
gence of such localization using the perturbation theory. 

In our work, we find that the extra link of the cycle 
indeed cause a different dynamical behavior compared to 
the dynamics on the cycle. The impact of the extra link 
is mainly determined by the largest eigenvalue and its 
corresponding eigenstate. This is similar to the other dy- 
namic processes on networks where the largest eigenvalue 
plays an important role in relevant dynamics. The addi- 
tional link causes a significant localization at the two end 
nodes of the link. We show such localization is related 
to the eigenstate of the largest eigenvalue, which differs 
from the ID regular cycles where the localization is only 
induced by the size of the system. This characteristic 
is reminiscent of the Anderson localization in condensed 
matter physics 



34| 



Anderson has shown that there is 
no quantum diffusion on disordered medium, an effect 
nowadays called (strong) localization, where the quan- 
tum mechanical transport through the lattice is prohib- 
ited by disordered structures js^,!!^. In our case, the ad- 
ditional link breaks the perfect order nature of the cycle 
and induces some structural disorder, this also leads to a 
strong localization resembling the Anderson localization. 
For Anderson localization all the eigenstates become lo- 
calized, whereas here it is only the one corresponding to 
the largest eigenvalue. We hope our findings provide a 
deeper understanding for the dynamics of quantum walks 
on networks and possible insights into the experimental 
implementation of continuous-time quantum walks on ir- 
regular networks (36| . 
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Appendix A: Definition of Chebyshev polynomials 

The Chebyshev polynomials of the first kind are defined by the recurrence relation [s^, [11] , 

To{x) = l,Ti{x) = x,2xT„{x) = r„_i(x) + T„+i(a;). (Al) 
The Chebyshev polynomials of the second kind are defined by the recurrence relation [s^, [13] , 

Uo{x) = l,Ui{x) = 2x,2xUn{x) = Un-i{x) + Un+i{x). (A2) 
The closed-form solutions of Eqs. (|A1I) and (IA2p are given by, 

Tn{x) = , (A3) 

— (n+1) _ 

Unix) - — , (A4) 



where z — x ~ y/x'^ — 1. For the case of large order n and z < —1 (jzj > 1), the above solutions can be approximated 

by, 

r„(x) « — , Unix) « -— ^ = -— (A5) 

2 \z — z\ z '- — z 

Using the closed- form solutions for the Chebyshev polynomials (See Eqs. (|A3|) and (|A4p ). we can prove the following 
identities, 

Un^iix) + U-n-iix)=0, (A6) 

2xUnix) = Un-lix) + Un+iix), (A7) 

Tnix) = Unix) - xUn-lix) = xUn-lix) - Un-2ix), (A8) 

Unix)Umix) - C/„_i(x)C/„_i(x) Un+niix), (A9) 

Unix)Umix) - Un+lix)Um-lix) = Un-raix), (AlO) 

U„{x)T„iix) + Um-iix)Tn+iix) = U„+mix), (All) 

Unix)Trnix) - i (x)T„+i (x) = Un-rnix), (A12) 

T,?(.t) - ix^ - l)U^_,ix) = 1, (A13) 

T„^ix)Unix) ^ -[Urn+nix) + Un-niix)], r„i(x)T„(x) = -[Tm+n{x) + r|„_„|(x)]. (A14) 

Appendix B: Calculation for the determinant equation 

For the sake of simplicity, we first simplify the four coefficients ci, C2, C3 and C4 using the Chebyshev polynomial 
identities. The first coefficient ci can be simplified as, 

ci = f/„_™(a:)[(2a; + l)U,n-2ix) - U^n-si^)] - U^^„i^iix)Um-2ix) 

^ {UM-mix)U„,^iix) - C/jv-m-i(a;)t^m-2(a;)} + U „^„,ix)U„,^2ix) (HZ]) uscd (Bl) 
= f/„_i(a;) -I- [/„_™(x)[/„_2(a;) dZSJ used 

The second coefficient can be simplified as, 

C2 = Ur,-niix)[Urn-4:ix) ~ (2x + l)Um-3ix) - 1] + C/jv-m- 1 (a;) J7to_3 (a;) - 1 

= [-Ur,-niix)Um-2ix) + C/jv-m- 1 (a^) (a^)] ^ ?7jv-m (a:) C/m-3 (a:) - U^-mix) - 1 (03) uscd (B2) 
= ~U^^2{x) - t/jv-m(a;)t/m-3(a;) - L/„_„(a;) - 1, dZS]) used 
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Analogously, C3 and C4 can be simplified as, 

C3 = UN-m-l{x)[{2x + l)U„i-2{x) - Um-3{x)] - U m -rn-2{x)Urn-2(x) + Um-2{x) + 1 

= {C/jv-m-l(a;)[/m_i(a;) - UM-m-2{x)Uyn^2[x)} + Ur,-rn-l{x)Urn-2ix) + Um-2{x) + 1 (B3) 
= [/jv-2(a;) + Ur,-m-l{x)Um-2{x) + Um-2{x) + 1, 

C4 = C/jv-m-i(a;)[C/„i-4(a;) - {2x + l)U,n-3{x) - 1] + C/jv-m-2(a;)J7m-3(a;) - Um-six) ~ {2x + 1) 

= {-U^-rn-l{x)Um-2{x) + C/jv-m-2 (a^) t/m-3 (2^) } - UM-m-l{x)Um-3{x) - C/«-m-l(a;) " Um-3{x) - (2x + 1) 

= -t/„_3(a;) - L/„_m_i(x)[/m_3(a;) - - [/„_3(x) - (2a; + 1) 

(B4) 

Substitute the simplified Ci, C2, C3 and C4 into C1C4 — C2C3 and expand the two products, we obtain, 

C1C4 - C2C3 =1 + UN~m{x) + U„i-2{x) - UN-lix) + {Ul^2ix) - Um-i{x)Uk^3{x)} 

+ {Um^2{x)U,^^„i{x) - Ur,-l{x)U^^,n^i{x)} + {U^-2{x)U„i-2{x) - i (x) C/„i_3 (x) } 

+ {Uf,-2{x)Urn-3{x) - U^-3{x)Um-2{x)}Uf,-m{x) + {U^-2{x)Urn-2{x) - i (x) t/™_3 (x) i (x) 

+ [U^-rn{x)Um-3{x) + UM-m-l{x)Um-2{x) - 2xUr,-rn{x)Um-2{x)] + 2[U^-2{x) - xUm-i{x)]. 

(B5) 

In the above equation, the terms in the big brackets { } can be simplified using (jAlOp . and the terms in the middle 
brackets [ ] can be simplified using (|A7|) and (|A8I) . Thus, 

C1C4 - C2C3 =2 + U^-„i{x) + 2U„i-2{x) - U^-i{x) + C/jv-m(a;) - C/jv_„i_i(a;)?7„_,„(x) + [/jv_m(x)C/jv_,„_i(x) ^^^^ 

+ [U^-rn{x)U„i-3{x) + UM-m-l{x)Um^2{x) - (C/jv_m-l(x) + t/„_,„+i (x)) C/,„_2 (x)] - 2Tjv(x). 



Applying identity (|A9p to the term in the middle bracket leads to, 

C1C4 - C2C3 = 2<^ 1 + t/„_„(x) + U^^2{x) - ?7jv_i(x) - T„(x) }> (B7) 



Using the Chebyshev polynomial identity (jA7[) in Appendix A, Eq. (1161) can be written as x^ = [Urn-i{x) + 
Um-2ix)]xm-i - [U,n-2{x) + [/m_3(x) + l]xm. According to Eq. (UHl), Xm-1 -C2Xm/ci and noting the simplified 
Cl, C2 in (IBip and (|B2p . thus we can write Xjv as. 



Xjv = ~{U,n-l{x) + Um-2{x))'^Xra ~ {Ura-2{x) + Ujn-3{x) + l)x™ 

t^«-2(a;) + C/jv-m(a;)[/m_3(x) + C/jv-m(x) + l . / ^ , rr / ^^ /rr ^ \ , rr ^ A , 1^ 
U^_^i^,) + U._,^{x)U^Mx) + - + + ^) 

t/iv-m(a;)[?7m-l(a;)C/„,_3(x) - ?7^t_2(x)] + [UN-2{x)Um-l{x) - UN~l{x)Um-2{x)] 

U,^-l{x) + U„-„i{x)Urn-2{x) 
_l_ [UN-2{x)Um-2{x) - C/w-i(x)[/„i-3(x)] + (x) i (x) + ?7„i-i(x) + U„i-2{x) - Um^i{x) 

Um-i{x) + Ur,-m{x)Um~2{x) 
m-l{x) + C/jv-m(a;) + C/jv-Tn(x)C/m_i(x) + Um-i{x) + t/„i_2(a;) - Uff-i{x) 
;7jv-l(x) + UM~7n{x)Um~2{x) 
_ t/w-m-l(a^) + UN-m{x)Um-l{x) + C/m-l(a;) + Um-2{x) - ^jv-l(x) ^ C2_^ 
Um-i{x) +UM-m{x)U^-2{x) " ~ Cl 

where the terms in the middle brackets ([ ]) are simplified using (jAlOp . 

Appendix C: Proof for the symmetry of the eigenstates 

In this section, we prove the symmetry of the components of the eigenstates, i.e., 

I I f \Xm+l-j\, if j G 

1 |x„+„+i_j|, if J e [m+ 1,7V] 



(B8) 



(Cl) 
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To prove the symmetric behavior of the eigenstates, we have the following theorem: If xi = ±Xm, is true, then all the 
other symmetric relations are true, i.e., 

If xi = ±Xm is true, then X2 = ±Xm-i, a^w = ia^m+i, •■• are also true. (C2) 



In the above theorem, the signs ± are uniformed. This theorem can be proved using Eq. (|12p and the Chebyshev 
identity. Set j = 2 and j = 3 in Eq. ([T2l) , we obtain, 

Xi = Urn-2{x)x„r^i - Um^3{x)Xm, (C3) 
X2 = Um-3{x)Xm-l ~ Um-4:{x)Xm- (C4) 

If Xi = ±Xm, then Xm = Um-2{x)xm-i/ (Urn-six) ± 1) (Eq. (jC3|) ). Substituting Xm into Eq. (jC4l) . we get the 
relationship between X2 and Xm-i, 

X2 TT I \ Urn-2{x)l}m-i{x) [Ul^_r^{x) - U,n-2{x)Uyn-A{x)] ±U,n-z{x) l±Urn-3{x) 

-Um-3[X) . . . = 77 . . . , = 77 , s , , = ±1, (^0) 



Xm-l U,n-?,{x)±l Um-3{x)±l Um-3{x)±l 

where the term in the middle brackets ([ ]) is simplified using (jA10[) . Thus we have proved X2 = ±a:„i_i. The 
remaining symmetric relations can also be proved in the same way (just use Eqs. (jl2l) and (|13p). The theorem 
suggests that to prove the symmetry of the eigenstates in Eq. (jCip , we only need to prove one symmetric component 
of the eigenstates {xi and x^) satisfies the Eq. (jCl|) . In the following, we will try to prove xi = ±x,„ using the 
determinant equation (|22l) . 

Now, we show that the determinant equation (j22l) can be written as a product of two factors using the Chebyshev 
polynomial technique. According to (jA12p . the term UN-m{x) can be written as UN~m{x) — U N~i{x)Tm-i{x) ~ 
U,ri-2{x)T^{x). Substituting Ufi^mix) into the determinant equation (j22p . we get 

{/„_i(x)T,„_i(x) = C/„_2(a:)[r„(a;) - 1] + r„(x) + C/„_i(x) - 1 (C6) 

Square both sides of the above equation and note that T^-^_^{x) — {x^ — l)U^_2{x) + 1 (Pell Equation in (jA13p ). we 
obtain, 

[{x' - l)Ul_,{x) - (r„(x) - l)2]C/,2„_2(^) ^2[T,{x) - WAx) + U.-i{x) - mn-2{x) 

+ [r„(x) - l][r„(a;) + 2U«-i{x) - 1] 

The factor in the middle bracket [ ] in the left hand side can be further simplified as, 

[{x^ -l)Ul_^{x)-{T^{x)-lf]^T^{x)~l~{T^{x)-lf = 2[T4x)-l]. dHH used (C8) 

Therefore, Eq. (jC7P becomes as, 

[T„(a:)-l]|2C/^_2(a^)-2[T„(a:) + t/„_i(a;)-l]f/„,_2(a:)-[r„(.T) + 2[/„_i(x)-l]l = [T^{x)~l]Q{N,m,x) ^ (C9) 



In the above equation, [Tn{x) — 1] = corresponds to the determinant equation for the cycles. Solving [Tn{x) — 1] = 
gives the eigenvalues for the regular cycle of N nodes, this suggests that some eigenvalues of the cycle are exactly the 
same as the eigenvalues in our model. Because we have squared the Eq. (jC6p , this doubles the number of solutions of 
the original determinant equation (l22t . [Tj^{x) — 1] = only gives about one half ([A^/2]) of the total solutions for 
Eq. (22), such eigenvalues belong to the original cycle and correspond to the eigenvalues whose first order corrections 
are (See Eq. ([291) 1 in the perturbation theory. In view of Fig. [TJ such eigenvalues are equivalent to the eigenvalues 
of a ID chain of size \N/2']. The other eigenvalues are determined by Q{N,m, x) — (the second factor in the big 
bracket) in Eq. ((C9l) . 

Next, we prove xi = using Eq. (|C9p . We use Eq. (fTS]) to write Xm-i = — C2Xm/ci, and substitute Xm-i into 
Eq. (IC3|) . 



Xl f.- , ,Ur,-m{x)Um-3{x) +U„-mix) +U,^-2{x) + 1 

— = <^ Um-2{x) rxTTr rvTr r~\ Um-six) 

[X Um-2[X 

_ [U,^-2{x)U„i-2{x) - U^-i{x)U„i^3{x)] + UN-m{x)Urn -2{x) + Urn-2ix) (CIO) 



;7jv_i(x) + U^-rn{x)Um-2ix) 
U^-niix) + U„^^{x)Ujn-2{x) + U„i-2{x) 



[/„_l(a;) + U,^-ra(x)Ujn-2{x) 



U^-l{x) + Ur,-mix)U,n-2{x) ^ (HIO)) USCd 
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If Tj^ix) -1 = 0, according to Eq. ((22)). U„-i{x) = U^-mix) + Uyn-2{x). Substitute it into Eq. (|C10|) . we get 
xxjxra = 1 (the denominator in Eq. (|C10p not equal to 0). If &{N,m,x) = (the second factor in the big bracket 
in Eq. (fCO]) ). according to Eq. ([22]), = [/jv_i(x) + Ti^{x) - Um^2{x) - 1, thus Ui^^rnix)Urn-2{x) can be 

simphfied as, 

(.t)C/„_2 (a;) = [U^-i{x)+T^{x) - U^-2{x) - l]U„,-2{x) = -]^[T^{x) +2U^-i{x) - 1], (Cll) 

where 8 (A^, TO, x) = has been used in the above calculation. Substituting L/N_„i(a;) + CAn-2 (a;) — C/jv-i(x)+r„(x) — 1 
(Eq. ([22])) and L/„_„(a;)L/™_2(cc) in ((CTT|) into Eq. ((CTO)) lead to, 

^ Un-i{x) + r„(x) - 1 - \[T,{x) + 2^^_i(x) - 1] ^ ^[r„(a:)-l] ^ _^ t fa:) ^ 1 fC12) 
Xm C/„-i(a;)- i[T„(x) + 2[/„_i(a;)-l] i[l-T„(x)] ' "^"^^ 

Thus we have proved xi = ixm, where the sign depends on the eigenvalues. For the largest eigenvalue, xi — —Xm 
holds because the largest eigenvalue satisfies Q{N, m,x) =0 in Eq. (|C9P {Tn{x) ^ 1). 

Now we prove that if m (or A'^ + to) is odd, the components Xj at the symmetry axis j = (m + 1)/2 of the eigenstates 
corresponding to eigenvalues satisfying Q{N,m,x) = 0, i.e., xj = (j = (to + l)/2) for the eigenstates whose 
eigenvalues satisfy 6(iV, to, x) = (See the second factor in the big bracket in Eq. (IC9p ). According to Eq. (IT^ . xj — 
U,n-i-j{x)Xjn-i{x)-U,n-2~]{x)x,n. Setting j = (m+ 1) /2 Icads to a;(„+i)/2 = Ui^,n-3)/2{x)x„i^i{x)-U(^ 

m— 5)/2 {,'^)'^m ■ 

Considering Xm-i — — C2^m/ci (See ci and C2 in Eq. (jBip and (jB2p ). we obtain, 



Ci=(™+l)/2 = S t^(m-3)/2(a;j ,(^\<TT I^^YTt 777\ U^„i-5)/2(x) ^X 



Uf,-rni x)Um-3{x) + UN-mjx) + Um^2 {x) + 1 
C/jv_l(x) + U,^-m{x)Um-2ix) 
UN~m{x)[Um~3{x)U(m-3)/2{x) - Um-2{x)U(^rn-5) /2{x)] + [UN-2{x)U(m-3)/2{x) - UN-l(x)U(,n-5)/2{x)] 

J7jv-i(.t) + U^-rnix)Urn-2ix) 

N m( ) (m 3)/2( ) (m 3)/2( (|A10[) cau bc uscd iu the middle brackcts 



J7jv-i(a;) + C/jv-m(x)C/,„_2(x) 

2{/jv-m(a;)C/(™_3)/2(a;) + Uj^^(^„i+l)/2{x) + U(rn-3)/2{x) _ KXr. 



J7jv-l(a;) + UN-m{x)Urn-2{x) Un-i{x) + UN^rn{x)Um-2{x) 

(C13) 

The numerator A in Eq. (IC13|) multiplied by r(m-i)/2(a;) and applying the identity (jA14l) . we obtain 
2r(„_i)/2(a;)J7(„_3)/2(a;) = Ujn-2{x) and r(m-i)/2(a:)?7iv-(m+i)/2(a:) = i[[/jv_i(a;) + C/jv-m(a;)]- Replacing the term 
Uff-m{x) = Un-i{x) + Tfi{x) — Um-2{x) — 1 (Scc Eq. (|22[) ) and after some appropriate algebraical calculations, we 
get, 

AT(„_i)/2(a;) = -U^^_2ix) + [T^ix) + U^-i{x) - l]Urn-2ix) + i[r„(x) +2;7„_i(x) - 1] = -^&iN,m,x). (C14) 

Compare Eq. (jC14|) and the second factor Q{N,m,x) in Eq. (jC9p . the above equation equals to if the eigenvalues 
satisfying Q{N,m,x) = 0. Thus A = and we have proved Xj = (j = (m + l)/2) for this case. Because the 
largest eigenvalue satisfies Q{N, m, x) — 0, the component x'j at the symmetry axis j = (to + l)/2 equals to for the 
corresponding eigenstate. The proof for j = {N + to + l)/2 {N + to is odd) is similar and we do not show here. 

Appendix D: Approximating the other eigenvalues using the determinant equation (|22p 

In this section, we approximate the other eigenvalues (except the largest eigenvalue) based on the determinant 
equation ([2^ . Eq. (|C6I) can be rewritten as, 

U„^-2{x)[T^{x) - 1] = U^^i{x)T„,^i{x) - T^{x) - U^-i{x) + 1 (Dl) 

Square both sides of the above equation and note that U^_2{x) — [T^_i{x) — l]/{x^ ~ 1) (Pell Equation in (jA13P ). 
we obtain, 

{{x^ - l)f/Li(x) - [T^{x) - lf}T^r^-lix) + 2\{x^ - l)ul_^{x) + (x^ - 1) (7„_i (x) [T^ (x) - 1]}t„_i(x) 

J ^ (D2) 

+ {{x^ - l)Ul_,{x) + 2{x^ - l)[/„_i(a;)[T„(x) - 1] + x^[T^{x) - = 
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Replacing the term (7^_i(a;) by ?7^_i(a;) = [Tl_^{x) - l]/{x'^ - 1) (See Pell EquationEHl leads to, 

[T^{x) - l]{2Tl_^{x) - 2[r„(x) + 1 + (x^ - l)U„-i{x)]Tm-i{x) 

^ (D3) 
+2{x^ - l)Ur,-i{x)+x^[T^{x) - l\+T^,{x) + 1 I = 

The first factor in the middle bracket [ ] in the left hand side T„{x^) — 1 = {Un-i{x^) = 0) corresponds to the 
eigenvalues for the regular cycle of nodes. Suppose the solutions of the second factor in the big bracket { } being 
equal to can be expanded as a; w a;'' — A (a;" is the solution of Tn{x^) — 1 = 0, A is a small value compared to a;°), 
the zero order approximation of Tm-i(a;), T„{x) and first order approximation of (x^ — 1)Un-i{x) are given by 

T„_i(x) «T„_i(xO), T„(x)«T„(x") = l (D4) 
- l)C/„_i(x) w (x°xO - l)[/„_i(xO) +xO[/„_i(x°)(x-xO) + 7VT„(xO)(x-x°) = -NA. (D5) 

Thus the second factor in the big bracket { } of Eq. (|D3|) can be approximated as, 

T^^^x") - (2 - 7VA)T„,_i(x") + (1 - NA) - 0, (D6) 

which leads to solution T,„_i(x°) = 1 — NA and the small deviation A = -^[1 — Tm-i(x^)]- Noting that r„i_i(x°) = 
cos(m — l)9n (the trigonometric function definition of the Chebyshev polynomials) [37j, |38|, our result obtained here 
accords with the first order correction in Eq. (I29p by the perturbation theory. 



M. Newman, A.-L. Barabasi, D.-J. Watts, The Struc- [20 
ture and Dynamics of Networks, Princeton University 
Press,2006. [21 
D.J. Watts, Small WoridsiThe Dynamics of Networks 
between Order and Randomness, Princeton University [22 
Press, 2003. [23 
S. Boccaletti, V. Latora, Y. Moreno, M. Chavez and D.- 
U. Hwang, Phys. Rep. 424, 175 (2006). [24 
A.-A. Juan and D.-J. Albert, New J. Phys. 9, 187 (2007). 
A. Alex, D.-J. Albert, K. Jurgen, M. Yamir and C.S. [25 
Zhou, Phys. Rep. 469, 93 (2008). [26 
P. S. Romualdo and V. Alessandro, Phys. Rev. Lett. 86, [27 
3200 (2001). 

J. G. Restrepo, E. Ott, and B. R. Hunt, Phys. Rev. Lett. [28 
100, 058701 (2008). [29 
M. E. J. Newman, Phys. Rev. E 74, 036104 (2006). 
J. Kempe, Contemp. Phys. 44, 307 (2002). [30 
V. Kendon, Math. Struct. Comp. Sci., 17,1169 (2006). [31 
O. Mulken and A. Blumen, Phys. Rep. 502, 37 (2011). [32 
N. Konno, Quantum walks, U. Franz and M. Schur- 
mann (Eds): Quantum Potential Theory, Lecture Notes 
in Mathematics, Vol. 1954, pp. 309-452, Springer, (2008). [33 
N. Shenvil, J. Kempe and K.-B. Whaley, Phys. Rev. A 
67, 052307 (2003). [34 
A. M. Childs and W. Dam, Rev. Mod. Phys., 82, 1 (2010). 
N. Ashwin and V. Ashvin, e-print [35 
|arXiv:quant-ph/0010117j 

N. Konno, Phys. Rev. E 72, 026113 (2005). [36 
O. Mulken and A. Blumen, Phys. Rev. E 71, 036128 
(2005). [37 
H. Krovi and T. A. Brun, Phys. Rev. A 73, 032341 (2006). 
O. Mulken and A. Blumen, Phys. Rev. E 71, 016101 
(2005). [38 



N. Konno, Inf. Dim. Anal. Quantum Probab. Rel. Topics 
9, 287-297 (2006). 

O. Mulken, V. Bierbaum, and A. Blumen, J. Chem. Phys. 
124, 124905 (2006). 

N. Konno, Int. J. Quantum Inf 4, 1023-1035 (2006). 

Y. Ide and N. Konno, Math. Struct. Comp. Sci 20, 1079- 

1090 (2010). 

J. Nagler, A. Levina and M. Timme, Nature Phys. 7, 265 
(2011). 

D. I. Tsomokos, Phys. Rev. A 83, 052315 (2011). 
X. P. Xu, Phys. Rev. E 77, 061127 (2008). 

M. D. Hatlee and J. J. Kozak, Phys. Rev. B 21, 1400 
(1980). 

T. C. Lubensky Phys. Rev. A 30, 2657 (1984). 

A. Blumen, J. Klafter, and G. Zumofen, Phys. Rev. B 

28, 6112 (1983). 

F. Jasch and A. Blumen, Phys. Rev. E 64, 066104 (2001). 

L. K. Gallos, Phys. Rev. E 70, 046116 (2004). 

O. Mulken, A. Blumen, T. Amthor, C. Giese, M. Reetz- 

Lamour, and M. Weidemller, Phys. Rev. Lett. 99, 090601 

(2007). 

E. Agliari, O. Muelken and A. Blumen, Int. J. Bifurcat. 
Chaos 20, 271 (2010). 

P. W. Anderson, Local moments and localized states. 
Rev. Mod. Phys., 50, 191-201 (1978). 
D.S. Wiersma, P. Bartolini, A. Lagendijk and R. Righini, 
Nature 390, 671-673 (1997). 

D. I. Tsomokos, S.Ashhab, and F.Nori, Phys. Rev. A 82, 
052311 (2010). 

J.C. Mason and D. C. Handscomb, Chebyshev Polynomi- 
als, Chapman and Hall/CRC; 1 edition (September 17, 
2002). 

T.J. Rivlin, Chebyshev Polynomials: From Approxima- 
tion Theory to Algebra and Number Theory, Wiley- 
Interscience; 2 edition (June 1990). 



